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I. INTRODUCTION 


The phenomenological predictions of snpersymmetry (SUSY) may be divided into three 
categories: (I) reflections of the supersymmetric Lagrangian in standard model phenomenol¬ 
ogy, including relations among the gauge coupling constants from SUSY grand unihcation 
and the presence of a heavy top quark and a light Higgs scalar; (II) the prediction of new 
particles with the correct spin and quantum number assignments to be superpartners of the 
standard model particles; and (III) well-defined quantitative relations among the couplings 
and masses of these new particles. While the predictions of (I) are of great interest, their 
verification is clearly no substitute for direct evidence. The discovery of a large number of 
particles in category (II) would be strong support for SUSY. On the other hand, the most 
compelling conhrmation of SUSY would likely be the precise verification of the relations of 
category (III). This would be especially true if, initially, only a small set of candidate SUSY 
partners are observed. 

Most discussions of supersymmetry at future high-energy colliders have concentrated 
single-mindedly on the question of particle searches. From one point of view, this is rea¬ 
sonable, because the existence of SUSY partners is unproven and this a prerequisite for any 
further analysis. On the other hand, the discovery of the first evidence for SUSY — or for 
any other theoretical extension of the standard model — will begin a program of detailed 
experimental investigation of the new sector of particles required by this extension. This 
investigation will need to be carried out with the same experimental tools that were used to 
make the original discovery. Thus, it is not only reasonable but also crucial, as we plan for 
the colliders of the next decade, to ask how any new physics that might be discovered can 
be examined in detail at these machines. 

Supersymmetry provides a particularly interesting subject for studies of the detailed 
analysis of physics beyond the standard model. SUSY models are weakly coupled, so their 
consequences can be worked out straightforwardly using perturbative computations. At the 
same time, SUSY models depend on a large number of unknown parameters, and different 
choices for these parameters yield qualitatively different realizations of possible new physics. 
Thus, the phenomenology of SUSY is quite complex. Eventually, if SUSY does give a correct 
model of Nature, the colliders of the next generation will be expected to determine the SUSY 
parameters, and their values will become clues that take us a step closer to a fundamental 
theory. We suggest that similar complexity should be found in any realistic extension of 
the standard model, and that similar investigations will be needed to understand the next, 
more fundamental, level. 

One consequence of the complexity of the parameter space of SUSY models is that it is not 
trivial to identify experimentally the specihc quantities which are related by supersymme¬ 
try. Faraggi, Hagelin, Kelley, and Nanopoulos |l|, Martin and Ramond @], and Kawamura, 
Murayama, and Yamaguchi have discussed in general terms the exploration of the spec¬ 
troscopy of supersymmetry partners, and the latter two groups have suggested particular 
mass relations which test supersymmetry independently of more detailed hypotheses. These 
tests are very ambitious, since they require mass measurements for the heaviest and most 
elusive particles of the superspectrum — the squarks, the heaviest partners of the Higgs and 
gauge bosons, and the sneutrino — at the 1% level. In these papers, very little attention was 
given to the question of how these experiments will be done. In this paper, we will present 


2 


some alternative tests of supersymmetry that involve only the lightest observable states of 
the superspectrum, and we will argue that these should be straightforward to carry out at 
colliders of the next generation. 

Our tests will exploit the advantages of the proposed Next Linear Collider (NLC), a 
linear e’''e“ collider with ^/s = 500 GeV and a design luminosity of 50fb“^/year [Q. This 
machine has already been shown to be a powerful tool for probing new physics In 

particular, previous work has shown that such a machine provides an excellent environment 
for measuring SUSY parameters under the assumption that newly discovered particles are 
sparticles . In this paper, we add to this body of work by showing 

how to test this assumption. Our analysis will take into account the relation of observable 
properties of the hnal state to the underlying reaction; as in the earlier NLC studies, we will 
be helped dramatically by the clean experimental environment expected at this machine. 
In addition, the expected availability of highly polarized electron beams should provide a 
powerful diagnostic tool. 

This study will be conducted in the context of the minimal supersymmetric standard 
model (MSSM). It is a reasonable expectation that charginos — the mixed superpartners of 
W bosons and charged Higgs bosons — will be among the lightest supersymmetric states, 
and that these will be accessible to the NLC. Thus, we concentrate here on tests of super- 
symmetry that involve the properties of charginos. The crucial problem we will face is that 
the mass eigenstates of charginos are in general a mixture of weak eigenstates, and their 
mixing pattern must be resolved before the quantitative implications of supersymmmetry 
become clear. To understand the experimental aspects of chargino reactions needed in this 
study, we have studied simulations of chargino production and decay using the parton-level 
Monte Carlo event generator of Feng and Strassler . 


The outline of this paper is as follows: In Sec. ^ we review the properties of charginos 
within the MSSM and state our assumptions. In Sec. 0 we divide the parameter space into 
characteristic regions. In Secs. |V| and |^, we present two different strategies for supersym¬ 
metry tests in two of these regions and analyze the experimental prospects for these tests 
in particular cases studes. In Sec. VI, we comment on other possible supersymmetry tests 
involving the properties of matter scalars and neutralinos. We present our conclusions in 
Sec. 


II. THE MSSM AND OUR ASSUMPTIONS 

Though our goal in the studies reported here is to test supersymmetry, we cannot be¬ 
gin without narrowing the phenomenological context. SUSY can, in principle, be realized in 
many ways. Here we assume that the observed particle content and qualitative phenomenol¬ 
ogy is that of the minimal supersymmetric extension of the standard model (MSSM), with 
conserved R-parity and therefore a stable lightest supersymmetric particle (LSP). This is 
the set of assumptions that is associated with the most commonly studied missing energy 
signatures for the discovery of candidate supersymmetric particles. R-parity conservation 
and the existence of only two Higgs doublets will be our two primary assumptions, and 
will be essential for much of the following analysis. We will also incorporate some minor 
additional restrictions for simplicity. In this section, we detail these assumptions and dehne 
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the basic set of parameters. A more detailed presentation of the MSSM can be found in 
many reviews |]A . 

The MSSM includes matter superhelds and two Higgs doublet superhelds Hi and H 2 , 
which give masses to the isospin — ^ and | particles, respectively. These two superhelds 
are coupled in the superpotential through the term —fj,eijH\H 2 , and the ratio of the two 
Higgs scalar vacuum expectation values is dehned to be tan/? = {H 2 )/{H'f). The MSSM 

mm 


also contains soft SUSY breaking terms 


which are parametrized by masses ruj for 


the scalar multiplets and masses Mi, M 2 , and M 3 for the U(l), SU(2), and SU(3) gauginos. 
In addition, there are cubic couplings (“A terms”) of Higgs scalars and sfermions. With 
the assumptions that we will make below, our study will be insensitive to the parameters 
entering through the A terms. 

The Higgsinos and electroweak gauginos of the MSSM mix to form two charginos and 
four neutralinos. In two-component spinor notation, the chargino mass eigenstates are = 
and = UijtlfJ, where = (—ilU^, H^) and, by convention, m-± < m^±. The 

matrices V and U diagonalize the mass terms 


{ip -I-h.c. 


( 1 ) 


where 


M^± = 


Mo 


a/2 M\y sin f3 


a/2 Mw cos /3 /i 


( 2 ) 


Ignoring some subtleties in this diagonalization having to do with negative mass values and 
the ordering of the eigenstates (see, for example, the first reference in |^), V and U are 


orthogonal matrices which can be parametrized by rotation angles 0+ and 0_. For (p± = 0, 
the chargino xf is pure gaugino, and for 0± = f, xf is pure Higgsino. The neutralino mass 


eigenstates are x? = 
mass terms 


where (ip^)'^ = {—iB,—iW^,H^,H 2 ) 


and N diagonalizes the 


-iryM^or+h.c. , 


( 3 ) 


where 


M^o 


( Ml 

0 

—Mz cos 
V Mz sin 


0 

I M 2 

iS [3 sin 9w Mz cos [3 cos 9w 
(3sm.9w —Mzsm(3cos9w 


—Mz cos [3 sin 9w Mz sin f3 sin 9w \ 
Mz cos [3 cos 9w —Mz sin [3 cos 9w 

0 —/i 

-p 0 / 


• ( 4 ) 


To reduce the large number of arbitrary parameters, we follow Ref. in introducing 
some additional assumptions. These assumptions are primarily phenomenologically moti¬ 
vated, and, where possible, we avoid assumptions based solely on grand unified theories 
(GUTs) and supergravity theories. As noted above, we assume R-parity conservation and 
the presence of a stable LSP, which we identify as the lightest neutralino Xi- In addition, 
we will ignore the intergenerational mixing in the quark and sfermion sectors, and we will 
assume that CP-violating phases in the SUSY parameters are negligible. We will also as¬ 
sume that one-loop effects do not introduce large and qualitatively new dependences on 


4 






















SUSY parameters. If these effects are large but may be absorbed by redefinitions of the 
tree level parameters, our analysis can be applied with only minor modihcations. The as¬ 
sumptions listed above will be in effect throughout this study. Additional conditions that 
are appropriate to the study of specihc processes and scenarios will be given below. 


III. THE PARAMETER SPACE OF CHARGINGS 


In many supersymmetric models, charginos are the lightest observable sparticles, and 
we now consider the possibilities for tests of SUSY from chargino production. As we are 
interested in what may be learned from the chargino signal, we will make, in this and the 
following two sections, the additional assumptions that gluinos, sfermions, and the Higgs 
scalars A°, and are beyond the kinematic reach of the NLC. Neutralino masses 
must be comparable to chargino masses, and below we will address the problem of removing 
neutralino backgrounds to the chargino signal. If a number of additional SUSY signals are 
available at NLC energies, their detection would be exciting in their own right, and would 
make possible the measurement of several sparticle masses. However, the procedure we 
outline below for measuring chargino couplings would not directly apply. Since we think it 
would be somewhat optimistic to expect a plethora of sparticles to be accessible at NLC 
energies, we have not explored this scenario further. 

The analysis of chargino pair production and decay is discussed in detail in Ref. JT^; 
here we will only summarize the most important qualitative features of this process. Using 
the picture of chargino production derived from this analysis, we will divide the parameter 
space into characteristic regions. In the following two sections, we will dehne and analyze 
tests of supersymmetry which rely on the particular characteristics of the chargino in each 
of these regions. 

Though the observables we will discuss involve only the chargino pair production cross 
section, the problems of experimental detection of the chargino signal necessarily bring in 
parameters of the chargino decay processes. We simplify our treatment of these processes in 
the following way: motivated by p —>■ ey and flavor changing neutral current constraints || 20 || , 
we assume that all left-handed sleptons of different generations are roughly degenerate (to 
within, say, 20 GeV) with mass m^, and the left-handed squarks of the hrst two generations 
are roughly degenerate with mass rrtq. In fact, chargino events are usually insensitive to all 
other sfermion masses. Decays through third generation squarks are suppressed because, 
for NLC energies, the mass difference m-± — m^o is almost always less than the top quark 
mass. For the remaining sfermions, the right-handed sfermion diagrams are suppressed by 
Higgs couplings ruf/Mw and are negligible. 

With these assumptions, there are only six parameters that enter the complete descrip¬ 
tion of chargino pair-production: /i, M 2 , tan/3, Mi, m^, and rriq. We do not assume gaugino 
mass unihcation, and so M 2 and Mi are unrelated. With an e~[ beam, XiXi production 
occurs through the s-channel Z and 7 diagrams and the t-channel z/g exchange diagram of 
Fig- 0, and so the left-handed differential cross section is governed by four parameters: 


dar 


dcos9 


CrC 


XiXi = 


) 


dar 


dcos9 


{fi, M 2 , tan/3, mi) 


( 5 ) 
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In the case of an beam, the z/g diagram is absent, and so the right-handed differential 
cross section is dependent on only the hrst three parameters: 

Charginos decay to the LSP either leptonically through W bosons or virtual sleptons, 

Xl ^ ^ x> , (7) 

or hadronically through W bosons or virtual squarks, 

xt Tq', W*) xlm', ( 8 ) 


and so all six parameters enter the decay process. The lighter chargino may also decay 
to LSPs through a virtual charged Higgs but this diagram is suppressed by Higgs 
couplings and is negligible for all but the most extreme choices of parameters. The heavier 
chargino may decay through complicated cascade decays. However, when xt production 
is kinematically accessible, the only information we will use about xt is its mass, which 
we will assume may be measured through threshold scans. The analysis will therefore be 
independent of xt branching fractions and other observables dependent on the details of 
the xt decay. 

The chargino masses m^± and m^± and the right-handed cross section an depend only 
on the parameters /r, M 2 , and tan/?, and these parameters may be used to dehne regions 
with qualitatively different behavior. To understand this, note hrst that, when M 2 3> |p| or 
|/i| S> M 2 , the following relations hold 
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m-± ~ min{|/i|,M 2 } and m-± ~ max{|/i|,M 2 } . 


( 9 ) 


These relations are in fact approximately valid in most of the available parameter space. 
The dependence of ctr on the parameters is more complicated. In Fig. ^ we plot contours of 
constant ctr for hxed tan (3 in the (p, M 2 ) plane. The dependence on tan (3 is fairly weak; we 
choose the representative value tan /? = 4 for illustration. Chargino production is inaccessible 
for y/s = 500 GeV in the hatched region, and the cross-hatched region is excluded by the 
current experimental mass limit m-± > 45 GeV |p2|,p^ . This leaves two bands, one on each 

Al ^ -• •- 


side of the p = 0 axis. At the top of each band, where M 2 3> |/i|, the chargino is Higgsino- 
like, xf ~ and we see that ctr is substantial. However as one moves into the region with 
M 2 ^ |/i|, (Tr quickly drops. This may be understood by noting that, because ^/s S> Mz, 
the 7 and Z production diagrams may be replaced to a good approximation by diagrams in 
which the U{1) and SU{2) gauge bosons B and are exchanged. However, the e/j couples 
only to B, and the couples only to W^. Thus, in the region with M 2 ^ |/i|, where the 
chargino is dominated by its wino component and xt ~ the cross section ur is highly 
suppressed. 

We are now in a position to dehne characteristic regions in the parameter space. These 
are shown for tan/? = 4 in the (/i, M 2 ) plane in Fig. The hatched and cross-hatched 
regions are as in Fig. In the remaining area, we dehne the following three regions, each 
of which includes a /i < 0 part and a corresponding /i > 0 part that is unlabeled: 
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Region 1: m^± < \fs- Here production is possible, and so both chargino masses 

can be measured. 

Region 2 (shaded): m-± + > ^/s, and < 10fb. 

Region 3 (shaded): m-± + m-± > -y/s, and aR ^ 50fb. 

These three regions almost completely £11 the region of parameter space in which chargino 
pair production is allowed at a 500 GeV e’''e“ collider, leaving only a small region in which 
the mixing is large and the chargino yf is just above threshold. In this study, we will ignore 
this small gap. In the two cases we will study in detail, we will assume m-± 172 GeV. For 

this value, the measurement of m-± constrains the parameters to lie on the dashed curves 
shown in Fig. |^. Then, if xt is not seen, ct/j < 10 fb or ctk > 70 fb for tan /5 > 4, and further, 
for 1 < tan /3 < 4, only small areas of the (p, M 2 ) plane he outside regions 1-3. For masses 
m-± nearer to threshold, the areas not covered by regions 1-3 are larger. However, this 
can be compensated by raising the collider center-of-mass energy, which increases the size 
of region 1. 

In region 3, if the ratio Mi/M 2 is fixed, x^f and become increasingly degenerate as M 2 
grows. Gharginos then decay to invisible LSPs and very soft jets and leptons. It is therefore 
difficult to choose a representative point in this region, as even the identification of the 
chargino signal can be difficult in some areas. More generally, if Mi and M 2 are unrelated 
(and, of course, independent of /i), m-± — m^o need not be small, even if the chargino 
is Higgsino-like. Although it may then be possible to verify SUSY relations in region 3, 
we will not consider this possibility further. However, we note that the MSSM makes a 
number of nontrivial predictions for region 3. Since xf ~ the u production diagram 
becomes negligible. The production forward-backward asymmetry is thus approximately 
zero. In addition, since the chargino is Higgsino-like, it decays predominantly through a 
virtual W, and so the ratio of hadronic to leptonic decays of the chargino should be equal 
to the corresponding ratio for W bosons. These characteristic features should distinguish a 
chargino candidate from new particles of other, non-supersymmetric origin. 


IV. A SUPERSYMMETRY TEST IN THE MIXED REGION 

We now study a representative point in region 1 in detail. The characteristic property 
of region 1 is that both chargino eigenstates can be produced, and so both masses are 
measurable. Thus, in this region, a promising approach will be to test the detailed form 
of the chargino mass matrix. In particular, notice that the matrix of Eq. (|^) contains, in 
addition to the new parameters M 2 , p, and tan/?, a dependence on the W mass. This is no 
accident. The off-diagonal matrix elements of Eq. (Q) result from the HWH vertex. This is 
related by supersymmetry to the HWdH vertex, which is related by gauge invariance to the 
term which gives mass to the W through the Higgs mechanism. Thus, verification that this 
parameter of Eq. (||) is indeed equal to Mw would be a quantitative test of supersymmetry. 
This test is formally independent of the neutralino sector and is therefore applicable to 
models with gauge singlets. 

We now investigate the extent to which we can realistically verify this correspondence at 
the NLG. In this example, and for the rest of this work, we will assume ^/s = 500 GeV. We 
will present results for integrated luminosities of 30 and 100 fb“^, corresponding roughly to 
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I to 2 years running at design luminosity. 

For our case study, we choose the underlying supersymmetry parameters to be 


(p, M 2 , tan/?, M 1 /M 2 , TUq) = (—195,210,4,0.5,400,700) . 
For these values, the MSSM gives 


( 10 ) 



m^± = 172 GeV 
m^o = 105 GeV 
= 255 GeV 


( 11 ) 


aR = 48 fb 
ai = 513 fb . 


For comparison, the QED production cross section is 397 fb. 

To investigate the expected sensitivity to the form for the chargino mass matrix, we 
generalize Eq. (^ to an arbitrary real 2x2 matrix, which we parametrize as 



( 12 ) 


Without SUSY, the ratio of off-diagonal elements need not be the ratio of vevs tan/? = 


{H 2 )/ (77i), and we have therefore replaced /? by /3^. As demanded by gauge invariance, we 
also replace Mz by = Mz{M^/Mw) in the neutralino mass matrix of Eq. (^. 0 

We will investigate to what extent the NLG experiments may conhrm the SUSY rela¬ 
tion M^ = Mw. More explicitly, we have extended the six-dimensional SUSY parameter 
space to a seven-dimensional parameter space, and we will investigate how well experiments 
may reduce the allowed region of this space to the supersymmetric subvolume in which 
M^ = Mw- Formally, this is a simple task. The four parameters entering Eq. (|T^) may be 
exchanged for the two masses and two mixing angles, (ni-±, nr-i, 0+, (/>_). By determining 
these four quantities from experiment, we can recover a constraint on M^. 

To determine the chargino masses and mixing angles from experiment, we will need to 
make assumptions about the decay properties of charginos. In our analysis, we will assume 
that these properties are those of a supersymmetric model at some point in parameter space, 
with the exception that the new chargino and neutralino mass matrices are used. Because 
we have not generalized the decay completely, this assumption is a compromise, but we feel, 
a reasonable one — it gives us a large but well-dehned space of possibilities to consider. 
In addition, we will see below, by explicitly scanning this space, that our results depend 
only weakly on the decay parameters. The main dependences are kinematic and would be 
expected in more general models of chargino decays. It is also worth noting that many of 
our assumptions may be checked a posteriori; for example, the assumption of a universal 


^The resulting neutralino mass matrix is not the most general allowed by gauge invariance. The 
fully general neutralino mass matrix will be considered briefly when neutralino events are considered 


in Sec. VI. 




left-handed slepton mass may be checked by observing the nniversality of leptonic branching 
fractions in chargino decay. 

The precision with which m-± and m^o can be determined was stndied by the JLC 
gronp 1^. Using a method that depends on kinematic argnments only, they fonnd that, for 
an integrated Inminosity of 20 fb~^, these masses conld be determined to approximately 2 
GeV, an nncertainty that is negligible for this stndy. The mass m-± may be determined 
by scanning near XiX^ threshold. Althongh a{e'^e~ —>• xtx^) is snppressed by abont an 
order of magnitnde from mixing angles, we will assnme that an energy scan will be able 
to determine m-± to a few GeV, and we will therefore also neglect this nncertainty in the 
following analysis. 

The crucial difficulty will be that of determining the two mixing angles. In principle, 
these can be extracted by measuring the right-polarized differential cross section for xt pair 
production, which is completely determined by the xt mass and the two mixing angles. The 
right-polarized cross section ct/j, though an order of magnitude smaller than ai, is still large 
enough to yield a sufficient number of events for precision studies. In particular, we will 
examine two quantities based on daji/d cos 9: the total cross section aji, and a truncated 
forward-backward asymmetry 

AX — '^Ri^ < cos^ < 0.755) — < cos 9 < 0) 

ctr(-1 < COS0 < 0.755) ’ ^ ^ 

where 9 is dehned as the angle between the beam and the positive chargino xt- (The 
motivation for this peculiar dehnition of will be given below.) With m^± known, the 
values of and determine the variables (0+, 0_) and may therefore bound M^. This 
strategy is appealing, because we have seemingly eliminated all dependence on three of the 
undetermined parameters of the theory: Mi, m^, and rriq. 

Unfortunately, the analysis is not independent of these three parameters when we con¬ 
sider what quantities are actually observable. Guts must be imposed to reduce standard 
model backgrounds. In this paper, we will rely on a standard set of cuts which have been pre¬ 
viously suggested to isolate the chargino pair production signal. These cuts select chargino 
events in which one chargino decays to an isolated hnal state lepton, and the other decays 
directly to hadrons. (Gharginos may also decay indirectly to hadrons through r leptons.) 
We will call such events “V mode events,” with the letter “V” chosen to suggest the typical 
2]+1 topology of these events. What is actually measured is not aR, but the Y mode partial 
cross section after cuts. 


rjay = 2rjBiBh(TR , (14) 

where r] is the efficiency of the cuts for Y events, Bh is the chargino branching ratio for 
direct hadronic decays, and Bi is the branching ratio for decays to a hnal-state lepton. 
These fractions both exclude decays to a r which subsequently decays hadronically. 

Since the charginos decay very quickly, with typical widths of 1-100 keV, the chargino 
direction and the asymmetry cannot be determined directly. We will measure A^ through 
its correlation to the forward-backward asymmetry of the hadronic system in Y events. 
In principle, the experimentally observable quantities AA^ and rjay depend on the decay 
distributions, and thus reintroduce dependence on the parameters Mi, and rriq. To 
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understand the extent of this problem, we have performed Monte Carlo simulations at a 
number of points in parameter space. These points have been chosen randomly, subject only 
to the constraints that they give values of m^o, and m-± consistent with those that 


X2 


would be measured in our case study. We will show below that, in the resulting subvolume 
of parameter space, the experimental observables turn out to be rather insensitive to Mi, 
mi, and mq, and therefore the virtues of our strategy in fact remain. 

To simulate chargino events, we used the parton level Monte Carlo event generator 
of Ref. 1^. This generator includes the spin correlations between production and decay 
processes. To simulate hadronization and detector effects, the hnal state partons were 
smeared with detector parameters as chosen in the JLC study H: 


—had 


40% 


and 


^lepton 


15% 


(15) 


where E is in GeV. 


The Y chargino events were selected by hrst using a system of cuts presented in Ref. [10 


These cuts are designed for charginos that decay through off-shell W bosons, and include 
the following: 

(a) I cos^jj < 0.9 for every hnal state parton, where 6i is the polar angle of parton i with 
respect to the e"*" beam axis. 

(b) El > 5 GeV, 9gi > 30°, that is, there must be an energetic e or /x with no hadronic 
activity within a cone of half angle 30°. 

(c) 20 GeV < invisible < a/s — 100 GeV. 

(d) ^acoplanarity ^ 150 . 

(e) ruhad < 68 GeV, i?had < \/s — 100 GeV, where ruhad and i?had are the mass and energy of 
the hadronic system. 

(f) [m^i, — Mil/1 > 10 GeV, where the u momentum is taken to be equal to the missing 
momentum. 

(g) —Qi cosOhad, Qe cos 9i < cos 41° = 0.755, where Qi is the charge of the isolated lepton, 
and 9i is as dehned in cut (a). 

These cuts isolate chargino events that have hadrons and an isolated lepton in the hnal state. 
We would like to isolate Y events, and we therefore need to eliminate events in which the 
hadronic system results from charginos decaying through r leptons. This may be done by 
imposing the additional requirement that the mass of the hadronic system mhad be greater 
than m^-. As was shown in Ref. [0, Y events very rarely have low ruhad at LEP II energies, 
and we have verihed that this is also true for NLG energies. We will therefore simply assume 
that this additional cut on ruhad cleanly isolates the Y mode events. 

Guts (c) and (d) are efficient for supersymmetric signals because of the large momentum 
and energy that are carried off by the unobserved massive LSPs. Guts (e)-(g) reduce the 
dominant standard model background, W pair production. In particular, cut (g) is designed 
to remove the large forward peak of WW events. Because the hadronic system’s polar angle 
distribution is truncated by cut (g), we choose as dehned in Eq. ([I^) , as the theoretical 
quantity with which we expect to be well-correlated. Since W pair production results 
primarily from annihilation, the use of these cuts in conjunction with a very highly 

right-polarized e~ beam results in a negligible background rate. The analysis of Ref. m 
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included ti events with a top quark mass of 150 GeV and found negligible background from 
this source. 

We caution the reader that the cuts (a)-(g) above have been designed to separate the 
chargino signal from standard model backgrounds, but have not been optimized to discrim¬ 
inate between xtxi production and other SUSY sources of Y events. In principle, these 
could include xfx 2 ^ind X 2 X 2 production, as well as the production of neutralino pairs 
XiXj- Ignoring effects of resolution smearing, the neutralino events will be backgrounds to 
Y events only when a heavy neutralino decays into a chargino and a W boson, which then 
decays leptonically to provide the single isolated lepton. While we have not simulated these 
events, we do not expect neutralinos to be a severe background because their production 
cross sections are generally small, and further, their decays to hxi and Zxi are usually fa¬ 
vored by phase space and therefore dominate. For the point that we are studying, the masses 
of the heavy neutralinos are m^o = 169 GeV, m^o =211 GeV, and m^o = 253 GeV. The 
decay X 4 is barely open, and the production of heavy chargino pairs is kinematically 

forbidden. Thus, xfx 2 production, with xf hF^Xi is the main SUSY contam¬ 

ination in the present case study. This background is restricted by phase space and mixing 
angles and can be eliminated entirely by running below the xfxJ production threshold. 

Throughout this study, we have assumed 100% beam polarization in our simulations. 
In the present case, however, because ai is an order of magnitude larger than the 
left-handed contamination of the right-handed beam could be substantial if the beam po¬ 
larization is not nearly 100%. If beam polarization near 100% is unobtainable, the signal 
may be determined by first measuring the signal to high accuracy, and then subtracting 
the left-handed contamination from the right-polarized e~ beam’s signal. For a beam polar¬ 
ization of 95%, these errors will not be large, and we have not included the statistical errors 
resulting from such a subtraction. It is clear, however, that highly polarized beams play a 
critical role in reducing such errors. 

We now determine the correlation of with through Monte Garlo simulations. 
A description of our method and the relevant formulae are contained in the appendix. We 
sample random points in the seven dimensional parameter space, with only the restriction 
that m^o, and m^± are each within 2 GeV of their values in Eq. ([Tl|) . For each set of 
parameters, we calculate from explicit analytical formulae and determine AA^ through 
Monte Garlo simulation. The results for 38 simulations are plotted in Fig. ^ A simple linear 
fit yields = 0.717A^ -|- 0.042 ± 0.036, where = 0.036 is the Icr deviation in 

for a fixed A\. The best £t is given by the solid line in Fig. and the la deviations are 
shown by the dashed lines. 

However, this quoted error overestimates the deviation from perfect correlation between 
A|. and because each point in Fig. ^ was computed from a hnite sample of Monte Garlo 
events and therefore contains a non-negligible statistical fluctuation. The average effective 
number of Monte Garlo events for the simulations was Nmc ~ 1400. Using the formulae 
contained in the appendix, we hnd that the Monte Garlo statistical error is AA^’^ = 0.026; 
when this is removed, the systematic error in assuming perfect correlation is found to be 
AA®^® _ 0.025. The correlation between A^ and AA^ is high — the chargino rest frames 
are slightly boosted, and the decay distributions are sufficiently similar for all sampled 
values of the underlying parameters that is highly insensitive to the decay process and 
is well-determined for a fixed A^. (If the beam energy is slightly reduced to run below 
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the xfx 2 threshold, the charginos will be less boosted. However, we do not expect the 
correlation between and to deteriorate much, since, even in the present case with 
■y/s = 500 GeV and only slightly boosted charginos, the correlation is high.) 

To determine the bounds that may be placed on experimentally, we must add the ex¬ 
perimental statistical error to For our representative point, a Monte Carlo simulation 

gives 


^had ^ -0.233 

T] = 35.5% (16) 

^exp = 6.0£/j , 

where Aexp is the number of Y events surviving the cuts, and Cr is the right-handed in¬ 
tegrated luminosity in fb~^. The total experimental uncertainties for two values of right- 
polarized integrated luminosity are found to be 

Cr = 30 (100) fb“^ ^ = -0.37 ± 0.107 (0.065) . (17) 

The efficiency rj also depends on the decay process. We determine rj by finding its 
range in the subvolume of parameter space in which the three masses and are within 
the experimental bounds of their underlying values. Each simulation gives a point in the 
(H^, rj) plane, and the distribution of points is plotted in Fig. A linear £t gives rj = 
—6.48A^ -|- 34.35 ± 1.07%, where = 1.07% is the la deviation in rj for a hxed A^. As 

in the previous hgure, the best linear £t is given by the solid line, and the dashed lines give 
the la deviations. We see that there is a dependence on — in cases in which chargino 
production is forward peaked, cut (g) lowers the efficiency. However, since we have already 
bounded A^ in the analysis above, we may use this measurement to restrict the range of rj. 
To determine the systematic error, we remove the Monte Carlo statistical error from Arj^Q. 
Following the analysis of the appendix, we hnd that Arj^ = 0.77% and Ar]®^® = 0.75%, 
and, including experimental statistical errors, we hnd 

Cr = 30 (100) fb"^ ^ = 8.0 (4.7) % . (18) 

ay 

To convert a measurement of ay into a measurement of aR, we must also take into 
account the uncertainty in the branching ratios Bi and B^. These again depend on the 
parameters of the chargino decay matrix elements and, in particular, on the masses 
and rriq. We have varied these masses to permit as a large a variation in aR as possible. 
However, the measurements of 'm^±, 'rn^±, A\, and ay constrain the allowed parameter 
ranges to regions where Xi ^md Xi have substantial Higgsino components. Recall that Bi 
and Bh take fixed values (equal to those for the W) in the Higgsino limit. These facts and 
the bounds m^, rriq > 250 GeV constrain the Y mode branching fraction to the region in 
which 29% < 2BiBh < 36%. Thus, the ay contours are rather insensitive to variations in 
the sfermion mass parameters. 

The measurements of A\ and aR constrain the plane to the shaded regions in 

Fig. The lightly (heavily) shaded region is the allowed region for Cr = 30 (100) fb“^. 
Contours of constant are also plotted in GeV, with the SUSY contour = Mw 
given by the dotted curves. The contours of constant aR that bound the allowed region run 
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roughly northwest to southeast; contours of constant run roughly southwest to northeast. 
The indicated boundaries correspond to Icr deviations in each quantity. 

Given the chargino masses of this case study, the theoretically possible range of is 


0<M^< 


m-± + m-± 

Xl X2 


218 GeV . 


(19) 


In the allowed region for Cr = 100 fb 

60 GeV <M^< 105 GeV . (20) 

The measurement of therefore, provides a quantitative conhrmation of SUSY. 

As an aside, we note that our analysis simultaneously bounds the parameters /r, M 2 , and 
tan/3^. In the heavily shaded region, the allowed ranges for these parameters are 

-204 GeV < /i < -183 GeV 

199 GeV < M 2 < 217 GeV (21) 

2.4 < tan/?^ . 


If one is led by the bounds on M^ (or other considerations) to view SUSY and the MSSM as 
conhrmed, one might then consider only the contour M^ = within the allowed region. 
One would also be led to identify tan with the ratio of Higgs scalar vevs, and so we will 
replace (3^ with f3. On the contour M^ = Mw, the bounds on the SUSY parameters are 
extremely strong: 

— 196 GeV < jji < —193 GeV 

208 GeV < M 2 < 211 GeV ( 22 ) 

3.9 < tan/3 < 4.1 . 


These bounds are so strong that it is likely that the uncertainties in chargino masses will 
be a significant source of uncertainty. (Recall that, while the uncertainties in chargino 
masses were included in the determination of systematic errors, the parameter bounds are 
determined from Fig. P, in which the chargino masses are hxed.) Nevertheless, it is clear 
that the discovery of both chargino mass eigenstates will allow one to place tight bounds 
on these three central SUSY parameters. In particular, the bound on tan/3 would be one 
of the most stringent and model-independent; the difficulty of determining tan/3 from the 
Higgs scalar sector is explained in Ref. m- Given the bounds of Eq. (p2|), other SUSY 
parameters may be restricted by additional measurements. For example, m^o may be used 
to determine Mi, and ai may be used to hnd m^. Such determinations may help lead us to 
an understanding of the SUSY breaking mechanism and other aspects of higher theories. 

We have now completed the case study for our chosen representative point. We conclude 
this section with comments concerning the power of this analysis for other points in region 
1. If one moves from the point given in Eq. ([T0| ) toward region 3, the results of the analysis 
become stronger for two reasons. First, an increases, and the experimental statistical errors 
decrease. Second, as a direct consequence of electroweak gauge invariance, such large values 
of aji can only be achieved for Higgsino-like even in the generalized (seven-dimensional) 
parameter space where one lets M^ vary. This implies that chargino decay is dominated by 
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the W diagram, and the sensitivity to the decay process parameters becomes even weaker 
than in onr case stndy. In particnlar, the systematic errors related to determining A\ and 
r] become smaller, and the branching ratios and take their W decay valnes. 

If one moves in the opposite direction toward region 2, the nnmber of right-polarized 
events deteriorates rapidly. In addition, xf may be gangino-like, and the branching fractions 
therefore depend more strongly on decay parameters, leading to a larger nncertainty in the 
determination of aR from ay. These problems can potentially be remedied by changing 
the analysis method. Since a highly right-polarized e~ beam leads to a very small level of 
backgronnd, it may be possible to nse a looser system of cnts, and to measnre the hadronic 
and leptonic branching fractions directly. The analysis in the gangino-like portion of region 
1 wonld then be limited only by statistics and systematic errors in the determination of A\ 
and rj, and the statistical nncertainties in the measnrements of the branching fractions for 
chargino decays. 

Finally, having considered variations of the Higgsino-gangino content of xf, one might 
consider variations orthogonal to these in the plane of Fig. |^, namely, variations in m-±. If 
xf is heavier, the chargino rest frame is less boosted relative to the lab frame. The decay 
process will then have a bigger effect on the correlation of AA^ with A\, and will 

increase. However, we have already considered a case with a fairly heavy xf: s-iid we see 
that the charginos need not be highly relativistic for to be small. In the opposite 

limit of lighter xf: the chargino rest frame is more boosted relative to the lab frame, decay 
effects become less important, and the resnlts of onr analysis can be expected to improve. 


V. A SUPERSYMMETRY TEST IN THE GAUGING REGION 

In the previons section, we considered the case in which both charginos were discovered, 
and fonnd that the SUSY constraint on the chargino mass matrix conid be verihed to 
fairly high precision. In this section, we examine region 2, in which only one chargino is 
seen and its prodnction cross section section from is small. Here we mnst rely on the 
chargino pair prodnction cross section from e^, which introdnces a strong dependence on 
mp from the second diagram in Fig. |I[ Fortnnately, there is an important compensating 
simplification: in this region, the charginos are very nearly pnre ganginos, and, in fact, it is 
a good approximation to neglect the deviations of cos0± from 1. In this limit, the conpling 
constant of the e^uxf vertex is related by snpersymmetry to the conpling constant 

g. Verihcation that this conpling constant is indeed eqnal to g wonld be a qnantitative test 
of snpersymmetry. 

For onr case stndy in region 2, we take the nnderlying snpersymmetry parameters to be 
(p, M 2 , tan (3, Ml/M 2 , m^-, m^) = (-500,170,4, 0.5,400, 700) . (23) 


For these valnes, the MSSM gives 
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( 24 ) 


m^o 

•4.1 

( 0 +, 0 -) 

(^R 

(Jl 


172 GeV 
86 GeV 
512 GeV 
( 1 . 2 °, 12 . 8 " 
0.15 fb 
612fb . 


For the point we have chosen (and for a signihcant part of region 2), the two-body chargino 
decay xf —>■ W'^Xi is open. The branching fractions Bi and are then fixed to their values 
in W decay, unless |/i| is very large, a possibility discussed at the end of this section. The 
case in which on-shell W decays are not allowed will also be discussed briefly at that point. 


To investigate the sensitivity of experiments to the value of the e^uxf coupling, we 
generalize this coupling from its SUSY value (/Vn to g^'Vu. We then test the SUSY relation 


= g. The differential cross section dai/dcosO is then a function of g 

0 and we can measure m-± 

Xi 


but because <;/)_ ~ 0 and we can measure m-±, we have only two unknowns. These may 

Xi 

be constrained with two quantities formed from dai/dcosO, which we choose to be ai and 


_ o'l{ 0 < cos6 < 0.707) — aL{—l < cos 6 < 0) 
“ (Ti(-1 < cos9 < 0707) 


(25) 


It is important to note that the parameters g^ and rriy enter dai/dcosO only through 
the v diagram amplitude, which has the form 


Ay ~ 




111 


t 


mt 


(26) 


Thus, for very large values of mjj, only the ratio g^/my can be determined. However, we will 
see that even for the rather large value of rriy that we have chosen, the two parameters g^ 
and rriy can be distinguished. In general, these parameters can be bounded independently 
when rriy is comparable to the collider center-of-mass energy (though still possibly above 
the pair-production threshold). 

We follow the procedure of the previous section, with the exception of using the cuts of 
Ref. p3, which are appropriate for charginos decaying through on-shell W bosons. These 


include the following: 

(a) > 5 GeV, 6qi > 60°. 

(b) j/t > 35 GeV. 

(c) ^acoplanarity ^ 150 . 

(d) — Mw\ > 10 GeV, where uisr is dehned to be the massless particle which, along 
with an initial state radiated photon in the Ez direction, makes up the missing momentum. 

(e) 6*sphericity < 45°, wliich we approximate in the Monte Garlo simulation by demanding 
—cos6 *had, QicosO^ < cos45° = 0.707. 

This system of cuts isolates chargino events containing hadrons and an isolated lepton. 
Again, the subset of these events that are Y mode events may be cleanly separated by 
demanding that ruhad be signihcantly larger than After these cuts, the WW background 
is reduced to roughly 25 fb for an e~[ beam, which is approximately the size of the signal after 
cuts. We will assume that the WW background is well-understood and may be subtracted 
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up to statistical fluctuations. As the WW background is strongly forward-peaked, we will 
also assume in computing statistical errors that it contributes completely to the set of events 
with COS0 > 0. The ti background, computed with = 140 GeV, is again negligible. In 
the gaugino region, the other SUSY signals do not provide a signihcant background to Y 
events, because the only kinematically accessible SUSY backgrounds are X 2 X 1 and X 2 X 2 ! 
with m^o ~ m^±. The neutralinos X 2 then decay to LSPs and an even number of leptons, 
and the number of events with one isolated lepton is highly suppressed. 

To determine the correlation between and A^, we perform Monte Carlo simulations 
at a number of randomly chosen points in the seven-dimensional parameter space (p, M 2 , 
tan /?, Ml/M 2 , mi, mq, g^), subject to the constraints that m-± and m^o are within 2 GeV of 
their measured values and ctr < 1 fb. Again the experimental observable is determined 
to be an excellent estimator of A^, with AA®^® = 0.034. A Monte Carlo simulation at our 
representative point gives 


^ Q_Q34 

T] = 11.9% (27) 

A^exp = 25.8Cl , 

where Cl is the left-handed integrated luminosity in fb~^. We now calculate the uncertain¬ 
ties in determining A^ and cxy using the equations found in the appendix, this time including 
also the errors arising from a substantial number of background events Yback ~ Yexp. We 
hnd that for two values of left-polarized integrated Inminosity, 

Cl = 30 (100) fb“^ ^ Al = 0.20 ± 0.067 (0.048) . (28) 

As in the previous case, the efficiency rj is fonnd to be highly constrained by the mea¬ 
surements of m^±, m^o, aL, and A^, and the resulting systematic error is = 0.55%. 

Including experimental statistical errors and those resulting from background subtraction, 
we hnd 

Cl = 30 (100) fb"^ ^ = 7.2 (5.6) % . (29) 

(Jy 

For Cl = 80 and 100 fb~^ these measurements constrain the allowed region of the 
(mp, gA plane to the shaded areas shown in Fig. |^. Because the charginos decay through 
on-shell W bosons, in contrast to the region 1 analysis, Bi and Bh are hxed at their values 
in W decay, and thus the contours for aL inferred from ay are independent of sfermion 
masses. For Cl = 100 fb~^, if mp < 250 GeV is excluded by the non-observation of any 
other threshold for heavy particle production, the allowed region is only the largest of the 
three shaded regions in Fig. |^. For this region, we hnd the constraint 

0.85^ <g^< 1.8g . (30) 

Such a result would be an important quantitative conhrmation of SUSY. 

Fig. I also illustrates a number of other interesting features. It is clear from Fig. that, 
without assuming SUSY, the analysis above has simnitaneously bounded the mass of a 
t-channel resonance, a useful result for future particle searches. If, on the other hand, we 
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assume the validity of SUSY, then we are restricted to the dotted line at /g = 1, and the 

measurement alone restricts rrij^. Alternatively, the ai measurement alone restricts 
to two different ranges, of which one can be immediately excluded. Finally, as expected from 
earlier comments, this analysis is signihcantly weakened if mp is large. For large m^, the 
contours of Fig. approach contours of constant g'^/my, and only the ratio g^/mc, can be 
determined. On the other hand, if it is possible to measure independently, for example, 
from e^z/yf production, then the bounds on g^/g can be signihcantly improved. 

In the example above, we have considered a point for which chargino decays through 
on-shell W bosons are allowed. This choice was motivated by two considerations. First, 
in region 1, we considered a point for which only off-shell W decays were possible, and 
appropriate cuts were used. Our choice in region 2 illustrates that tests of SUSY are also 
possible when cuts appropriate to on-shell W decays must be used. Second, the scenario in 
which on-shell W decays are possible becomes more and more typical as the chargino mass 
rises, and the analysis presented is thus generalizable to higher chargino masses and beam 
energies. It is easy, however, to hnd points in region 2 where the chargino cannot decay to 
an on-shell W. For example, if one assumes the GUT relation M 2 = 2Mi, on-shell W decays 
are excluded for m^± ^ 160 GeV. In this case, we must use the cuts presented in Sec. B 
In addition, chargino decays through virtual sfermions are not negligible, and one must 
consider the dependences of the branching ratios on sfermion masses. Such dependences 
will introduce systematic errors that may considerably weaken our results. However, as in 
the case of the gaugino portion of region 1, if these branching ratios can be measured, the 
systematic errors in their determination may be greatly reduced. In contrast to the region 
1 case, the e/ beam, with its accompanying WW background, must be used. However, 
because WW events do not usually have ]/t without isolated leptons, they are likely to be 
a small background to purely hadronic chargino events. Although further study is required, 
it again seems probable that the Y mode branching fraction can be measured directly, and, 
with these modihcations, the previous analysis may be applied to region 2 scenarios in which 
only off-shell W decays are allowed. 

It is also true that in the very far gaugino region with |/r| S> M 2 , where yf ~ W^ 
and ~ B, the W decay diagram is suppressed by mixing angles, and, even when decays 
through on-shell W bosons are kinematically allowed, virtual sfermion diagrams may be 
important. This requires that xf and Xi be very nearly pure gauginos, however, and this 
occurs only for |/i| > 1 TeV, a condition that is disfavored by hne-tuning constraints. 


VI. SFERMIONS AND NEUTRALINOS 

Up to this point, we have considered only precision SUSY tests from studies of the 
properties of charginos. Other sparticles may be produced at NLG energies, however, and 
we now examine the possibility of testing SUSY through the properties of sfermions and 
neutralinos. The discussion will be limited to brief remarks and, in contrast to the previous 
sections, no attempt will be made to perform detailed studies. 

We hrst investigate the possibility of identifying a few newly-discovered scalars as 
sfermions. We are most interested in the scenario in which these scalars provide the hrst 
opportunity for precision tests of SUSY, and we therefore consider the case in which these 
scalars are lighter than charginos. In contrast to the previous sections, we will not impose 
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any constraints on intergenerational slepton and sqnark mass degeneracies. However, if the 
problem is considered in fnll generality, it is complicated by many arbitrary parameters 
associated with sfermion intergenerational mixing. Simply to make the problem tractable, 
we will assnme that intergenerational mixing is absent. We will also assnme that left-right 
mixings may be neglected, with the nnderstanding that the discnssion that follows may not 
be applicable to the sfermions of the third generation. Probes of the left-right mixing of 
scalar tans have recently been discnssed by Nojiri . 

With these assnmptions, the properties of these sfermions are completely specihed by 
their qnantnm nnmbers and their masses. The only category (III) tests involving sfermion 
properties are therefore verihcations of mass relations. Given the assnmptions above, the 
masses of sfermions are 


^UL = 
2 

ml = 

dh 

= 

2 

rrij = 

m?^ = 
= 

= 


rriQ + 101,1 + ~ I sin^ 6w) cos 2j3 

rriQ + m1 + M|(—| | sin^ 6w) cos 2f3 

m| + ml + M|(| sin^ 6w] 




'd + Mzi 


i sin^ 6w 


cos 2/3 
cos 2(3 


(31) 


m| + ml + M|(—I -f sin^ 9w) cos 2/3 

m| -I- |M| cos 2/3 

mi, + m‘( + Mli— sin^ 9u/) cos 2/3 . 


where m^, mp, mp, mp, and mp are soft SUSY breaking scalar masses. Similar relations 
hold for second and third generation sfermions. With additional relations from grand nni- 
hcation, there are a nnmber of relations among these scalar masses [Q. However, we will 
continne to eschew assnmptions that are not phenomenologically motivated. Withont GUT 
assnmptions, the right-handed masses are nnrelated to the other masses, and the left-handed 
masses are related only by 


m|^ - m|^ = cos 2(3 

“ - Mil,, cos 2(3 , 


(32) 


where we have omitted the small fermion mass terms. For tan/3 > 1, these mass differences 
are positive, bnt we will consider all possible valnes of tan (3 below. The relations of Eq. (|32D 
are qnantitative predictions of SUSY that we may try to test. 

Unfortnnately, if the newly discovered scalars are sleptons, it will be impossible to test 
these relations, becanse the snentrinos will decay invisibly throngh vl —> vx\. We are as- 
snming that charginos are heavier than snentrinos, and so the decay v —*• eixl is exclnded. 
Also, even if tan (3 <1 and m^^ > mg^, the experimental lower bonnd mg^ > 45 GeV p2|j23[ 
implies mp^ — rrie^ < 47GeV, and the decay vl —>• clW^ is also greatly snppressed. The 
masses of sleptons are therefore highly nnlikely to provide the first category (HI) verihca¬ 
tions of SUSY. Of conrse, if snentrinos are heavier than charginos, precise verihcations of 
slepton mass relations conld be nsed to snpplement measnrements of chargino properties. It 
shonid also be noted that other properties of slepton events may provide additional preci¬ 
sion measnrements in the gangino region. It may be possible, for example, to measnre some 
nentralino properties throngh the 3-channel exchange diagrams for charged slepton pair 
prodnction 0. 

On the other hand, if the scalars are sqnarks, both left-handed species will decay visibly. 
A previons stndy of sqnark mass determination fonnd that at the NLG, in most regions 
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of parameter space, squark masses can be measured to approximately 2 GeV with an in¬ 
tegrated luminosity of 10fb“^, even in scenarios with cascade decays |T^. This study also 
found that left-handed squarks can be effectively separated from right-handed squarks us¬ 
ing beam polarization. It may be difficult to properly assign flavors to the different squark 
mass thresholds, however, especially if these thresholds are not well-separated. Let us hrst 
suppose that the masses of only two left-handed squarks are determined. To verify SUSY 
quantitatively, one must assume that the squarks are in the same generation, and must also 
independently determine tan/? from the Higgs scalar sector. This is by no means always 
possible, and most likely requires, for example, that m^o < 300 GeV so that a heavy Higgs 
boson is kinematically accessible 1^^. Even if all of these measurements can be made, the 
precision of the test is not high. For example, if rriq > 200 GeV, the mass difference is 
< 15 GeV, and so in the best case scenario where tan/3 is determined exactly, 
the squark mass relation can be verihed to approximately 20%. If it is not possible to mea¬ 
sure tan (3 from the Higgs boson sector, a precision test of squark mass relations is only 
possible if one measures four left-handed squark masses. One can then check that there 
exists some flavor assignment consistent with 


"^1 - ~ ’ ( 33 ) 

where |A(m^)| < M^. 

The possibility of making the hrst quantitative tests of SUSY from sfermion properties 
is therefore not very promising. In the case of sleptons, the prospects are bleak, while in the 
case of squarks, even after assuming that intergenerational mixing is absent, precision tests 
are complicated by difficulties in havor determination and rely on many MSSM scalars being 
kinematically accessible. However, the sfermion sector provides a number of opportunities 
for disproving the MSSM and SUSY. For example, if sneutrino decay is observed, one of 
our assumptions must be invalid. Also, the relations of Eq. (|32D are valid not just for the 
MSSM, but are extremely general predictions of SUSY. If they are found to be violated, not 
only will the MSSM be excluded, but almost all supersymmetric models will be strongly 
disfavored. On the other hand, of SUSY is favored by experiment, measurements of the 
squark and slepton masses will give important information about the havor dependence of 
the SUSY breaking mechanism. 

Neutralinos are natural candidates for precision SUSY tests, because, with the assump¬ 
tion that the lightest neutralino is the LSP, all sparticle event observables depend, at 
least formally, on the parameters that determine neutralino properties. In addition, neu¬ 
tralinos are light in many models, and, in fact, throughout parameter space, if charginos are 
produced, X 1 X 2 production is kinematically possible. 

One might hope to follow the procedure in Sec. Q by generalizing the neutralino mass 
matrix. If we relax SUSY, the most general form of Eq. (|^) consistent with gauge invariance 
is 



/ 

V 


Ml 0 
0 M 2 


M ^ 

0 —fi 
-/i 0 


where M is an arbitrary 2x2 matrix that may be parametrized as 


( 34 ) 
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( 35 ) 


_ (—Mz COS (3^ sinsin 6 *^ \ 

y Mz cos j3^ cos 9^ —sin/?^ cos^^^ j 

There are then seven parameters that enter neutralino events, and one must try to check 
the SUSY relations 9^ = 9w and = 1. A general analysis is likely to 

be complicated. One possible simplihcation would be to consider a less than fully general 
neutralino mass matrix by setting, for example, = 1. On the other hand, one might wish 
to assume the standard SUSY neutralino mass matrix, generalize the neutralino-fermion- 
sfermion coupling to , and check that = g- However, even this analysis is more 
complicated than the chargino case, because the SUSY neutralino mass matrix contains 
an additional parameter. In addition, an important caveat to all analyses based on the 
neutralino mass matrix is that such analyses rely on the absence of gauge singlets, and are 
therefore more model-dependent than the chargino analyses of previous sections. 

Without detailed study, it is not possible to dismiss the possibility that precision studies 
of sfermion and neutralino properties may be useful for testing SUSY. However, even from 
the brief comments presented above, it is clear that the sfermion and neutralino sectors 
are signihcantly less promising than the chargino sector. Category (HI) tests from chargino 
properties are likely to be the least model dependent and may be the hrst strong quantitative 
tests even if some other sparticles are lighter than charginos. 


VII. CONCLUSIONS 

Softly broken supersymmetric theories are like spontaneously broken gauge theories in 
that the relationships between dimensionless couplings implied by the symmetry continue to 
be preserved, while the corresponding relationships between the masses of various particles 
can be badly violated. It is this feature which provides the best opportunity for quantitative 
tests of supersymmetry. In this study we have examined the possibilities for testing various 
SUSY relations in a number of scenarios. These studies have been conducted in the experi¬ 
mental setting provided by a linear e’''e“ collider with polarizable beams, and results have 
been presented for = 500 GeV and integrated luminosities of 30 and 100 fb“^. 

In the scenario in which charginos are the hrst sparticles to be discovered, we have ana¬ 
lyzed two representative cases. In the hrst, we probed the form of the chargino mass matrix, 
and in the second, we tested the Xiff coupling. In both examples, we found that the test 
led to rather strong quantitative conhrmations of the MSSM and SUSY. As a by-product, 
interesting bounds on some SUSY parameters were also obtained. The availability of po¬ 
larizable beams was found to play a vital role, allowing us to dehne characteristic regions, 
ehectively eliminate dependences on certain SUSY parameters, and remove background. 
Our analysis was performed using a parton level Monte Carlo event generator and did not 
incorporate possible contamination of chargino pair events from other SUSY processes. Of 
course, a more detailed analysis that includes the simultaneous production of all possible 
SUSY events together with a more realistic simulation is needed before dehnitive conclusions 
about precision SUSY tests may be drawn. 

The prospects for obtaining the hrst quantitative tests of SUSY from sfermion and neu¬ 
tralino properties were also considered. Sleptons were found to be poor candidates for such 
tests because of the difficulty in detecting sneutrinos, and precision tests from squarks were 
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found to rely on the discovery of at least four squarks or two squarks and, most likely, two 
Higgs bosons. The analysis of neutralino properties is complicated by its dependence on 
a large number of parameters. Whether these complications may be overcome in certain 
scenarios remains to be seen in further studies. However, while falsihcation of sfermion mass 
relations is the least model-dependent disproof of SUSY, it is likely that the chargino sector 
is the simplest and most powerful for verifying the quantitative predictions of SUSY. 

We have not considered the possibilities for quantitative SUSY tests at other colliders, 
nor have we examined the additional constraints that come with the adoption of GUT and 
supergravity assumptions. Even with fairly weak assumptions, however, we have found 
that, if sparticles are produced at future e+e“ colliders, measurements of their properties 
may allow us to quantitatively verify SUSY, a valuable first step in the exploration of the 
full structure of supersymmetric theories. 
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APPENDIX: UNCERTAINTY ANALYSIS 

In this study we use the truncated forward-backward asymmetry = Af, where i = L 
or R, and the Y mode partial cross section ay to constrain parameter space. These the¬ 
oretical quantities are found through their correlations to experimental observables. The 
uncertainties in determining A^ and ay therefore receive contributions from two sources: 
systematic errors, that is, uncertainties arising from the lack of perfect correlation between 
the theoretical quantities and the experimental observables, and experimental statistical er¬ 
rors. In this appendix we collect the formulae used to estimate the systematic and statistical 
errors. 

Systematic errors are determined by performing Monte Carlo simulations at a number 
of points in parameter space. The truncated forward-backward asymmetry of chargino 
production before cuts, is determined through its correlation to the forward- 

backward asymmetry of the hadronic system’s direction after cuts. The theoretical quantity 
A^ depends only on parameters that enter the production process, while depends on 
both production and decay, and on cuts and detector effects. The systematic uncertainty 
in A^ is therefore determined by the sensitivity of to the decay process, cuts, and 
detector effects, and this sensitivity is measured through simulations. For each of iVpts 
points in parameter space, A^ is determined from exact analytical expressions, and A*^®^ 
is found from a Monte Carlo simulation. A linear fit to the resulting distribution in the 
(A^, plane is parametrized by 

A*^®^ = aA^ + b± AA'j^^c , (Al) 
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where is the la uncertainty in for a hxed A^. The total Monte Carlo uncertainty 

includes both the systematic error and fluctuations from hnite Monte Carlo statistics. 
The contribution from Monte Carlo statistical fluctuations is 


where 
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(A3) 


is the Monte Carlo statistical uncertainty in A^*^^ for simulation i, and Nua is the effective 
number of events in simulation i. The systematic error of the distribution is then 


= TAisyMTifif. (a4) 

To the systematic error must be added the experimental statistical error. This error is 
given by 
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where A^^d jg ^he forward-backward asymmetry for our case study, and Aexp (A^back) is the 
number of signal (background) events that pass all cuts and is proportional to the integrated 
luminosity. (Here we have assumed that the background is well-understood and may be 
subtracted up to statistical uncertainties. We also assume that all background events are 
in the forward hemisphere, a good approximation for the dominant background, W pair 
production.) We estimate the total experimental uncertainty in for a given A^ to be 

AA*°; = {AA^y^f + [AAf-^y . (A6) 

What we actually measure is however. We therefore are more interested in the exper¬ 
imental uncertainty in A^ for a hxed which is 


AA>^ = |a|-'AA“ , (A7) 

where a is the slope of the linear £t in Eq. ( |A1|) . 

The efficiency of the cuts r] is found simply by its correlation to previous measurements. 
To determine the uncertainty in rj, we reduce the parameter space to the region in which 
the previous measurements have their appropriate values and determine the variation of rj 
within this subspace. We determine rj for each of the simulations and obtain a distribution 
of points in the (A^, r/) plane. The best linear £t to this distribution is 


r, = a'A^ + h'± , 


(AS) 
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where A? 7 ^q is the Icr error in rj for a hxed A^. To hnd the systematic error, we must again 
remove the fluctuations that arise solely from hnite Monte Carlo statistics. The Monte Carlo 
statistical error is 



where C is the integrated luminosity. For the purposes of this study, /S.C/C is negligible. The 
uncertainty in the number of Y events passing the cuts is AiVexp = y^Aexp + Aback, where 
Aexp (Aback) is the number of Y mode (background) events passing the cuts, respectively, 
and we have again assumed that the background is well-understood and may be subtracted 
up to statistical uncertainties. 
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FIGURES 

FIG. 1. The diagrams contributing to chargino production at e'’'e“ colliders. The t-channel 
diagram is absent for beams. 

FIG. 2. Contours of constant (in fb) for fixed tan/3 = 4 in the (/r,M 2 ) plane. Chargino 
production is inaccessible for ^fs = 500 GeV in the hatched region, and the cross-hatched region 
is excluded by the current experimental mass limit > 45 GeV. The cross section oji quickly 
drops to zero in the \fi\ ^ M 2 regions. 

FIG. 3. The three characteristic regions for fixed tan/3 = 4 in the (/r,M 2 ) plane, as defined 
in the text. (The corresponding /r > 0 parts of these regions are unlabeled.) The hatched and 
cross-hatched regions are as in Fig. 0 The dashed curve is the contour m-± = 172 GeV. 


FIG. 4. The correlation of and for 38 points in the seven-dimensional parameter space 
(/r, M 2 , tan/3^, Ml, m^, rUg, M^). These points have been picked randomly, subject only to the 
constraints that 'm^±, , and are within 2 GeV of their underlying values in the case study. 

The linear best fit is given by the solid line, and the la deviations are given by the dashed lines. 


FIG. 5. The correlation of rj and for the 38 points in the seven-dimensional parameter space 
(/r, M 2 , tan /3^, rriq, M^), selected as in Fig. Q. The linear best fit is given by the solid line, 

and the la deviations are given by the dashed lines. 


FIG. 6. The allowed region of the (</>+, 4>-) plane from measurements of A^ and ay- The lightly 
(heavily) shaded region is allowed for Cji = 30(100) fb“^. Gontours of constant M^ are plotted 
in GeV. On the dotted contours, the SUSY relation M^ = Mw holds. 

FIG. 7. Allowed regions (shaded) of the plane for Cl = (a) 30 fb“^ and (b) 100 fb“^. 

Solid (dashed) curves are contours of constant aL (^x) that bound the allowed regions. On the 
dotted lines, the SUSY relation = g is satisfied. 
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